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There has been a significant progress in the last decade in deriving precisely, 
from lattice calculations, the static interaction potential between two heavy quarks. 
In the so-called quenched approximation (only valence, no sea quarks to start with) 
a funnel potential containing a linear plus a color-coulomb term is well established™ 
Spin and velocity dependent corrections to this form have been obtained as well^^. 
When sea quarks are incorporated (unquenched approximation) the long-distance 
behavior may change dramatically. As a matter of fact it has been shown in QCD 
at finite temperature 1^ and also in SU(2) Yang-Mills theory El that the potential 
saturates, i.e., it gets a constant value from a certain distance. Physically the sat- 
uration of the potential is related to screening: light qq pairs are created out of 
the vacuum between the two heavy-quark sources giving rise to a screening of their 
color charges. The effect coming out at short distances is that the running of the 
QCD coupling slows down with the distance. This translates in an effective coulomb 
strength bigger than in the quenched case. At long distances (bigger than the satu- 
ration distance) string breaking may take place: it becomes energetically favorable 
the recombination of the light q and q of the pair with the two heavy sources to 
form two static-light mesons. 

Although saturation has not been proven in QCD D one can wonder whether 
there is any direct indication of it from real QCD and/or phenomenology. We shall 
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argue that this question has a positive answer. Then we shaU construct a potential 
model incorporating screening for the description of hadrons. Its application to the 
non-relativistic bottomonium system will allow us to study the phenomenological 
effect of screening and to predict a lower bound for the saturation distance. 

From real QCD the Schwinger-Dyson equation allows to extract the running 
coupling, as, behavior in terms of an effective gluon mass, 

Mg(Q2) This be- 
havior is connected to confinement. Actually when this coupling, particularized 
for Mg(Cp — > 0) ~ ^QCD, is implemented into the one-gluon-exchange diagram a 
confining linear potential comes out. More generally a certain connection can be un- 
veiled by assuming a Yukawa parametrization of the confining potential as proposed 
long time ago from lattice simulations!^ 

Vconfir) = err = or . (1) 

\ fir J fi fir 

By realizing that the linear behavior is recovered for fi, the inverse of the satu- 
ration distance, going to 0, the following connection has been proposed 1^ 

KQ^) = Aqcd - MgiQ^) . (2) 

Then plausible values for the saturation distance, for heavy meson scales, go 
from /i-i(Q2 = (5 GeV)^) = 1.25 fm to fi'^iQ^ = (1 GeV)^) = 2.15 fm. 

From phenomenology the finite number of experimentally detected hadronic 
resonances could well be the reflection of the saturation property of the poten- 
tial. Indeed from a calculation of the non-strange baryon spectra with a Yukawa- 
screened interaction a perfect one-to-one correspondence between quark model 
bound states and experimental resonances can be obtained El For heavy mesons 
the same parametrization allows for an adequate description of bb and to a lesser 
extent of cc. However this smooth parametrization does not seem to correspond to 
more recent lattice indications of a rather abrupt saturation transition ^ . To check 
the effect of such a fast deconfinement we shall analyze its consequences in the 
simplest non-relativistic system, bottomonium, by using the potential model: 

Vir) = { ^<^- , (3) 

1^ arbr ~ e/rbr r > rbr 

where tr is a confining parameter whose value can be guessed from the value of the 
force at intermediate distances ^'^ (cr ^ 800 — 900 MeV-fm~^), e is the unquenched 
Coulomb strength (e = 106 MeV-fm (a^ — 0.4)), and rbr is the breaking or sat- 
uration distance for which we expect 1 fm < r^r < 2 fm at the {Q'^)bb scale. In 
order to be more precise we fit the bb spectrum. By choosing the b quark mass 
altogether with a so as to reproduce the ground state mass and spectral energy 
separations (mf, = 4820.5 MeV, a — 800 MeV-fm~^), rbr is fitted by assuming that 
the energy breaking or saturation threshold, Ebr, is just above the highest energy 
resonance known, say Ebr 11050 MeV. In this manner a lower bound value for rbr 
{rbr = 1-76 fm) is in fact obtained. The results for the complete spectrum (up to d 
states) are shown in Table 1 where the calculated masses are assigned to spin-triplet 



2, 2008 0:54 WSPC/INSTRUCTION FILE p-gonzalez 



bb description with a screened potential 3 

states. Since differences in masses with data are at most of 30 MeV an unambiguous 
identification of states is possible though for the ns states with n > 4 some room 
for a slight mixing with d states is left. 



Table 1. bb spin-triplet bound state masses and properties. 





Mass (MeV) 


Exp. 




(r2)i/2 (fm) 


Is 


9460 


9460.30 ± 0.26 


0.10 


0.21 


2s 


10022 


10023.26 ±0.31 


0.08 


0.51 


Id 


10171 


10162.2 ± 1.6 


0.02 


0.56 


3s 


10346 


10355.2 ±0.5 


0.09 


0.76 


2d 


10444 




0.04 


0.80 


is 


10605 


10580.0 ±3.5 


0.10 


0.97 


3d 


10679 




0.06 


1.00 


5s 


10829 


10865 ± 8 


0.11 


1.12 




10889 




0.07 


1.23 


es 


11030 


11019 ± 8 


0.11 


1.43 


Ip 


9935 


9900.1 ± 0.5 


0.03 


0.41 


2p 


10267 


10260.0 ±0.5 


0.05 


0.67 


3p 


10531 




0.07 


0.89 


4p 


10761 




0.08 


1.09 


5p 


10958 




0.08 


1.43 



This is confirmed by the calculation of leptonic widths, Table 2 (for details of 
the calculation see Ref. [TJ, showing a good agreement with experiment (see the 
fourth column where the effect of radiative and relativistic corrections is partially 
eliminated for ns states by taking their ratio to the Is case). 

Table 2. Leptonic widths F^+g- for bb in keV. 



rfi^_{i-^) (r,+,-),,p rfl^_/rf^^_(is) [r,+,-/r,+,-(is)] 



Is 


1.01 




1.32 ±0.05 


1 




1 


2s 


0.35 




0.520 ±0.032 


0.35 




0.41 ± 0.07 


Id 


2.2 X 10" 


-4 




2.2 X 10- 


-4 




3s 


0.25 




seen 


0.25 




seen 


2d 


2.6 X 10" 


-4 




2.6 X 10" 


-4 




4s 


0.22 




0.248 ± 0.031 


0.22 




0.19 ± 0.03 


3d 


5.7 X 10" 


-4 




5.7 X 10" 


-4 




5s 


0.18 




0.31 ±0.07 


0.18 




0.24 ±0.06 


4d 


6.7 X 10" 


-4 




6.7 X 10" 


-4 




6s 


0.14 




0.130 ±0.030 


0.14 




0.10 ±0.03 



Electromagnetic El transitions have been also evaluated, Table 3, in excellent 
agreement with data. 

Finally the spin singlet-triplet splitting has been calculated perturbatively. Our 
result M[T(ls)] — M [rih{ls)] ~ 185.7 MeV agrees with the experimental value 
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Table 3. El decay widths for bb in keV. 



Transition 




Texp 


T(2s) - 




1.62 


1.7±0.5 


T(2s) - 


-7X61 (IP) 


2.55 


3.0 ± 0.8 


T(2s) - 


-7X62 (IP) 


2.51 


3.1 ±0.8 


T(3s) - 
T(3s) - 
T(3s) - 


-7Xf.„(2P) 
-7Xi,i(2P) 
-7X62(2^') 


1.77 
2.88 
3.14 


1.4 ±0.4 
3.0 ±0.6 
3.0 ±0.6 



M [T(ls)] - M [77f,(ls)] = 160 ± 40 MeV within the experimental errors. 

We can then conclude that our naive non-relativistic quark model provides with 
an excellent description of bottomonium, much better than the corresponding to 
a non-screened funnel potential and at the level of the best descriptions obtained 
with much more ellaborated models. These results seem to support the view that a 
rather abrupt saturation may take place in real QCD. Moreover we predict that at 
the momentum scale corresponding to bb the saturation distance should be bigger 
than 1.75 fm. From this large value one does not expect saturation can be easily seen 
in current lattice simulations with the usual Wilson loops techniques. Meanwhile 
it could be worthwhile to examine at the phcnomenological level the effects of an 
abrupt saturation in other meson and baryon systems. 
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